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Abstract 

In quantum theory particles are represented as wave packets. Shock 
wave analysis of quantum equations of motion shows that wave function 
representation in general and wave packet description in particular con¬ 
tains discontinuities due to a non-zero quantum force. The quantum force 
causes wave packet dispersion which results in the intersection of char¬ 
acteristic curves developing a shock discontinuity. Since quantum force 
vanishes for localized quantum density waves [1], it is thus established 
that localized quantum density waves form the only class of wave func¬ 
tion representation of particles in quantum theory without shock wave 
discontinuities. 


1 Introduction 

Wave-particle duality is a central aspect of quantum processes. It is well known 
that the Schrodinger equation describes such a duality in the form of wave 
packets solutions, particularly with the Gaussian wave packets [2-5]. Wave 
packet solutions of the Schrodinger equation possess an interesting properly 
that even in the absence of an external potential they exhibit dispersion. Such 
dispersion occurs at very short time scales, causing interference for electron 
waves. 

In quantum potential formalism, wave packet dispersion results from the ten¬ 
dency of quantum trajectories to accelerate away from each other. Quantum tra¬ 
jectories for a free Gaussian wave packet, given by x{t) = u^t+XQ^Jl + (/it/2mcrg)^, 
indicate that components of the wave packet starting off at close by yet different 
initial positions and initial speeds intersect after a short time, exhibiting a shock 
formation. Intersection of characteristics lines is a typical aspect of nonlinear 
systems. However, in quantum dynamics of wave packets, quantum trajectories 
are assumed to be non-crossing. In the following exact analysis of quantum wave 
packet dynamics we drop this assumption, and the complete nonlinear system 
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is analyzed using general theory of Riemann invariants. This leads to the in¬ 
teresting result that quantum wave packets develop shock wave discontinuities 
immediately after their formation. 

Beginning with the next Section, we investigate the existence of shock wave 
phenomenon in quantum hydrodynamic formulation of the Schrodinger wave 
equation. It is shown that for a general wave function solution of the Schrodinger 
equation a non-zero quantum force causes characteristics to intersect, hence 
generates shock wave discontinuities in a quantized system. Such discontinu¬ 
ities have a travelling wave form, and correspond to particle motion in the free 
particle case. We take the example of a Gaussian wave packet to calculate the 
position and time of formation of quantum shocks in Section 3, whereas Section 
4 gives a summary of the main conclusions of the work and its relation to some 
recent experiments on electron waves. 


2 Shock Wave Analysis 


In the quantum potential approach, the general form of the wave function can be 
written as V'(r, t) = i?(r, t) exp iS'(r, t)//i. Then the time-dependent Schrodinger 
equation, with an external potential: ihdtp{r,t)/dt = —+ 
V{r)'ip(r,t) gives, after separating real and imaginary parts, the following equa¬ 
tions 

+ V-(p(r,t)v(r,t)) = 0, (I) 

+ (v(r,t) • V) v(r,t) = -—V{V{r)+mQ{r,t)). (2) 

ot m 

where Q{r,t) = —(fi^/2m^i?(r, t))V^i?(r,t), and p(r,t) = Equations 

(1) and (2) are the basic equations of quantum dynamics [6] in a fixed Eulerian 
frame, with respect to which the relative velocity of an element is v(r,<). 

We write the equation (I) and (2) as a single matrix equation, and keep to 
the one dimensional case only. Then equations (I) and (2) give. 


r 1 

0 ■ 


o 

I 



u p 
Q p ^ 



Ux 


( 3 ) 


where u is the component of the velocity v along the x-direction. Here pt = 
dp/dt and Qp = dQ/dp, etc. 

The eigenvalues and eigenvectors for system (3) can be calculated from the 
characteristic equation 


det „ ^ = 0, 

Qp u - \ 

( 4 ) 

which gives 


{u - Xf - pQp = 0. 

( 5 ) 

Thus the two eigenvalues Ai _2 are 


Ai,2 =u± y^pQp. 

(6) 
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According to Riemann theory of shock waves [7], these eigenvalues give the 
characteristic speed for families of characteristics C+ and C-. Thus we have 


dx 


C+ : \ — u + a /pQp — Ai, 

dx 


C_ : A = — = « - = A2, 

and the corresponding eigenvectors are given by 


■ Ml ' 


^y^pIQp 

. _ 


1 


( 7 ) 

( 8 ) 

( 9 ) 


In general, equations for characteristic lines for the system can be written 
as 

X±(t,to) = a;(fo) + {t + to), (10) 

where {dx/dt)^^^^ = \r{t = to)- For shocks to develop, these characteristic (10) 
must intersect at some common point in space. This occurs if the slope of each 
characteristic increases with the initial time to- In view of equations (7) and (8) 
this is the case provided Qp ^ 0, that is, if the quantum force dQ/dx does not 
vanish. Another way to state shock condition is to expand function X±{t,to) 
for 6to < < 1 in the neighborhood of each characteristics as 

X±(t,to + Sto) X±(t,to) -\ -- ^-Sto, (11) 

Oto 


then any two neighboring characteristics X±{t,to) and X±{t,to + Sto) intersect 
provided 


dX±(t,to) 

dto 


= 0 . 


( 12 ) 


This is the shock condition for the system (3), which we shall use in the following 
to calculate the time of shock formation. 

Having obtained the necessary (and sufficient) shock conditions (12), we can 
now explicitly determine the (travelling) shock wave solutions for the system 
(3). This is done in the Appendix, using Riemann invariants. These calculations 
show that the shock waves develop if the quantum force is non-zero, thus the 
phenomenon is of purely quantum nature. Since for a free Gaussian wave packet 
the quantum force is non-zero in general, this implies the existence of quantum 
shocks in the case of free Gaussian wave packets. 


3 Quantum Shocks for the Case of a Gaussian 
Wave Packet 

Gaussian wave packet solution to the Schrodinger equation for a free particle is 
given by the wave function 

'!/'(a;, t) = -exp (ik{x — uot/2) — {x — uot)^) /4.scro, (13) 

(27rs)'’' 


3 



where uq is the uniform constant speed of the wave packet, and the measure of 
the spread | s \= a = \/(Jq(1 + {ht/2ma'lY). Correspondingly, the amplitude 
and the phase functions are given by 

t) = -— exp -{x - (14) 

(27rcr^) ' 

S(x, t) = —^ tan~^( ^ y t) + mu(x — uot/2) + (x — uot^ ^ „ t, (15) 

2 2mcjQ amaQa^ 

and the quantum potential is, 


Q{x,t) 


dm^cr^ 2cr2 


(16) 


Using u = (dS/dx) /m, we obtain from equation (16) the speed of a wave packet 
element 


u = uo + 




dm^CTncr^ 


{x - Uot). 


(17) 


hence by integration, the position of a wave packet component is 


x{t) = Uot + Xo 




^2man 


(18) 


Here xq and uo denote the initial position and velocity of the wave components, 
respectively. 

Since by the above analysis the characteristics along which each wave packet 
component travels must intersect, we determine the location of the quantum 
shock in this case. Using equations (7) and (8), the equation of characteristic is 
given by 

-^±(t,to) = x(to) + (u± y'pQf}] {t + to). (19) 

\ / t —to 

Then from the shock condition (12), we have 


d{u± \/pQp)^_^ / , _\ 

x{to)-{ - — - —(t + to) - (u± ^JpQp] =0. (20) 

ato \ ' / t=to 

Substituting for Q, u, and x from equations (16), (17) and (18) respectively, 
and then taking to = 0 we have, after some simplification. 


ts 


—8to^CTq 

K^xo 



and therefore, from equations (7) and (8), 


Xs 


{uo ± 



( 21 ) 


( 22 ) 


Equations (21) and (22) give the time and position for the quantum shock 
associated with the Gaussian wave packet (13). 
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4 Conclusions 


The analysis presented in this paper shows that if in any region of space the 
quantum force tends to increase, a shock-like situation will develop. This must 
be so, since slope of the characteristics (19) then increases, causing character¬ 
istic curves to intersect. This was explicitly shown for the case of Gaussian 
wave packets, where the quantum force is equal to h^{x — uot)/4ma{t)'^; which 
although decreasing first, then increases, and eventually attains a constant limit 
as a; —>■ oo and t —>■ oo (Fig. 1). This indicates that the quantum force causes 
the wave packet to burst rather than spread smoothly. 

If the relative velocity of the Gaussian wave packet and the lab frame 
coincides, quantum shock occurs at time tg = \J 4m^CTQ^/fi^Xg, and position 
Xg = jxQ, travelling with the speed Xg/tg. Equation (22) shows that 

this speed differs from the classical formula by a constant H/ due to 

wave packet dispersion. 

A similar analysis can be applied to the case of Airy beams [8], and to the 
recently observed leviton structures [9], where limits on electron interferom¬ 
etry has been reduced to the attosecond scale [10,11]. For quantum density 
soliton waves, representing particle-like localization, such discontinuities do not 
form, since quantum force in this case is identically zero. This result also has 
implications for the problem of equivalence principle in quantum theory [12]. 

APPENDIX: Riemann Invariants and the Shock Wave Solution 

To obtain the Riemann invariants for the system (3), we take a linear com¬ 
bination of the equations (3) with coefficients being the components of the 
eigenvectors (9), this gives: 


{pt + UPx + pUx) ± \JP/Qp {Ut + QpPx + UUx) = 0. (Al) 

Since {dp/dp) / p/Q^ = {du/dp), it follows that 


du Qp dp 

dp y p du \j Qp 

(A2) 

Then substituting for Ux = UpPx = ^^Qp/ppx, and px = PuUx 
into equation (Al) we have 

= i\/ p/QpUx, 

[pt -b ± yJpQp'j Px'j ± \J p/Qp {ut + (u± y'pQp'j Ux^ 

= 0, (A3) 

and using equations (7) and (8): 


d{u± F{p)) = 0, where F{p) = J \/Qp/pdp- 

(A4) 
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Integrating we have the two constants (Riemann invariants), A and B along the 
characteristics: 


u + F{p) = A(^), along Fi with parameter ^ = x + \t, (A5) 

u — F{p) = B{ri), along r 2 with parameter rj = x — Xt. (A6) 

Thus eliminating u and F(p), we have the solution for u{x,t) and p{x,t)' 

u{x, t) = A(x + At) + B{x — At), (A7) 

F{p) = A{x + At) — B{x — At). (A8) 


which can be easily verified as a (traveling) shock wave solution to the system 

(3). We notice in the above analysis that if the quantum force is zero, then the 
shock speeds Ai and A 2 equal the speed u. 
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Eigure Caption: 

Eig. 1: Quantum force for the Gaussian wave packet (13), with h, m, ao 
unity, and uq = 10 units. 
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